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I.  INTRODUC  ION 

This  report  is  the  final  one  of  a  series  describing  our  investigation  of 
the  effects  of  surface  roughness  on  emissivity  and,  in  particular,  of  the  re- 
suits  of  our  calculation  of  the  microwave  brightness  temperatures  of  the  sea. 

Prior  to  this  study  the  on’y  theory  of  rough  surface  emissivity  was  that  of 

1  2 
Stogryn  who  used,  in  the  reake  representation  for  the  emissivity,  a  single- 

scatter  geometrical-optics  approximation  for  the  surface-scattered  electro¬ 
magnetic  fields.  To  the  resulting  equation  for  the  thermal  energy  emitted  by 
the  sea  surface  was  added  the  microwave  energy  originating  in  the  atmosphere 
and  scattered  (also  in  the  single  scatter  approximation)  by  the  sea  surface. 

The  total  energy  leaving  the  surface  (described  in  terms  of  an  "apparent"  or 
"brightness"  temperature)  was  evaluated  by  Stogryn  for  a  variety  of  wind  speeds 
and  observation  angles.  Stogryn's  calculations  provided  the  first  insight  into 

the  dependence  of  radiometric  measurements  on  s»a  state,  polarization,  and  ob- 

3 

servation  angle.  However,  it  js  readily  shown  that  Stogryn's  theory  gives 
erroneous  results  for  surfaces  with  appreciable  roughness.  Use  of  the  simple 
single-scatter  approximation  in  the  theory  has  been  shown  to  lead  to  a  non¬ 
physical  creation  or  loss  of  energy  with  the  result  that  calculated  emlssivities 
and  sea  brightness  temperatures  may  differ  significantly  from  the  correct  values. 

The  source  of  the  inconsistency  in  the  Stogryn  method  was  identified  and 

a  correct  geometrical  optics  formulation  was  developed,  first  for  a  cylindrical 

roughness  model  of  the  sea  surface  and  subsequently  for  a  general  two- 

c 

dimensional  roughness  model.  It  was  shown  in  this  series  of  reports  by  the 
authors  that:  (1)  the  single  scatter  approximation  is  inadequate  and  it  is 
essential  to  Include  both  surface  shadowing  and  (at  least)  double  scatter  by  the 
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randomly  rough  surface;  (2)  the  corrected  theory,  in  marked  contrast  to  the 
Stogryn  theory,  conserves  energy  to  a  high  degree  of  approximation;  (3)  the  geo¬ 
metrical  optics  theory,  including  shadowing  and  double  scatter,  provides  rigorous 
upper  and  lower  bounds  to  the  effect  of  surface  roughness  on  the  emissivity  and 
brightness  temperature;  (k)  application  of  the  theory  to  the  sea  showed  that, 
for  wavelengths  such  that  geometrical  optics  is  appropriate,  the  effect  of 
changing  sea  state  on  radiometric  measurements  can  be  calculated  to  a  satisfac¬ 
torily  high  degree  of  accuracy.  The  model  used  accommodates  changing  water 
temperature  and  salinity,  varying  atmospheric  conditions,  any  sea  state,  all 
observation  angles,  and  both  polarizations.  An  exhaustive  series  of  computer 
calculations  were  carried  out  to  establish  the  dependence  of  measured  sea 
brightness  temperature  on  all  the  possible  variables.**  The  effects  of  foam  and 
spray  on  total  brightness  temperature  were  excluded  from  the  model  because  of 
the  lack  of  reliable  experimental  or  theoretical  data.  Approximate  semi- 
empirical  models  for  the  effects  of  foam  and  spray  can  be  readily  incorporated, 
however,  as  was  illustrated  in  the  case  of  the  earlier  cylindrical  roughness 
mode  1 . 


Comparison  of  the  calculated  sea  brightness  temperatures  with  experimental 

6 

values  showed  good  agreement  for  wavelengths  in  the  neighborhood  of  1  cm,  while 
for  longer  wavelengths  discrepancies  appeared.  This  is  to  be  expected,  for  the 
strict  validity  of  geometrical  optic  requires  the  surface  to  have  negligible 
curvature  over  distances  of  the  order  of  a  wavelength;  with  increasing  wave¬ 
length  the  surface  will  be  Increasingly  structured  over  wavelength-sized 
dimensions  and  deviations  from  the  purely  geometrical  optics  predictions  would 
result.  The  qualitative  effect  of  surface  height  variations  within  a  wave¬ 
length  (which  we  will  refer  to  as  "small-scale  structure")  can  be  predicted/ 


<  * 
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The  effect  on  a  radiation  field  of  surface  structure  of  a  given  size  decreases 
with  Increasing  wavelength;  Indeed,  for  sufficiently  long  wavelengths  roughness 
effects  on  scattered  or  emitted  radiation  must  disappear  entirely.  Thus,  for 
longer  wavelengths,  an  increasing  fraction  of  the  rough  surface  will  radiate  as 
though  it  were  smooth  and  the  total  roughness  effect  must  diminish  (assuming 
there  is  no  strong  periodic  component  in  the  height  spectrum  which  would  pro¬ 
duce  resonance  effects).  This  general  diminution  of  the  roughness  dependence 
of  the  emissivity  is  enhanced  by  the  fact  that  the  roughness  parameter  for  com¬ 
ponents  of  surface  structure  with  scales  larger  than  the  radiation  wavelength 
is  the  mean-square  slope,  and  for  the  case  of  the  sea  surface  the  mean-square 
slope  associated  with  such  large  scale  components  is  known  to  decrease  with 
increasing  wavelength.  This  expected  general  decrease  in  sensitivity  of  the 

sea  brightness  temperature  to  roughness,  with  decreasing  frequency,  was  borne 

8  9 

out  by  the  experiments  of  Hoi  linger  *  for  frequencies  between  1.4  and  19.4  gHz 
and  wind  speeds  of  up  to  14  m/s. 

There  are  some  additional  theoretical  and  experimental  arguments  which 
have  been  advanced  regarding  the  dependence  of  the  sea  brightness  temperature 
on  sea  state  for  longer  microwave  wavelengths.  The  possibility  exists  that 
the  geometrical  optics  model  will  continue  to  apply,  but  with  the  roughness 
parameter  —  the  rms  slope  of  the  total  surface  —  being  replaced  by  a  smaller 
"effective"  rms  slope  which  is  wavelength  dependent.  This  possibility,  as 
well  as  the  form  of  the  effective  roughness  parameter,  was  suggested  by  Lynch  ^ 
who  formulated  a  theory  of  rough  surface  scattering  for  wavelengths  sufficiently 
long  that  significant  surface  curvature  may  exist  within  dimensions  of  the 
order  of  a  wavelength  He  found  that,  except  for  near-grazing  angles,  the 
scattered  intensity  is  still  described  by  the  geometrical  optics  theory  but 
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with  the  actual  rms  surface  slope  replaced  by  a  reduced  rms  slope  whose  value 
depended  on  both  wavelength  and  observation  angle;  the  same  behavior  would  carry 

i 

over  to  the  emission  problem.  Similar  conclusions  were  suggested  by  Hoi  linger 
who  showed  that  his  data  for  the  change  in  brightness  temperature  with  wind  speed, 
at  a  fixed  observation  angle,  could  be  fit  by  the  geometrical  optics  theory  pro¬ 
vided  one  assumed  an  effective  rms  slope  for  the  surface  whose  value  decreased 
with  increasing  wavelength. 

However,  to  relate  sea  brightness  temperatures  to  widely  varying  sea  sur¬ 
face,  atmospheric,  and  measurement  conditions  with  sufficient  accuracy  to  permit 
the  use  of  a  passive  microwave  system  as  a  viable  tool  for  remote  measurement  of 
sea  surface  conditions,  a  more  precise  theory  or  model  of  rough  surface  emis¬ 
sivity  is  required.  The  geometrical  optics  theory,  developed  in  our  earlier 
reports,  is  evidently  adequate  for  frequencies  of  about  20  gHz  and  higher,  but 
for  lower  frequencies  a  more  general  theory  is  needed.  The  required  generali¬ 
zation,  which  includes  the  effects  of  both  large  and  small-scale  surface 
structure  on  sea  brightness  temperature,  is  described  in  this  report.  In  Section 
II  we  discuss  our  general  approach  to  the  problem,  the  definition  of  the  com¬ 
posite  surface  model,  and  the  form  of  the  wave  height  spectrum.  In  Section  III 
the  theory  of  the  emissivity  of  an  anisotropic,  slightly  rough  random  surface 
is  applied  to  a  modeling  of  the  small-scale  sea-surface  structure  and  some 
numerical  results  are  described.  The  equations  used  are  a  new  re.  jsentation 
for  the  emissivity,  correct  through  second  order  in  the  roughness  parameter, 
and  are  derived  in  the  Appendix.  Section  IV  contains  a  summary  of  the  numerical 
results  using  the  complete  composite  surface  model  of  the  sea.  The  mode1  is 
shown  to  conserve  energy  at  least  as  well  as  the  geometrical  optics  model,  and 
the  relation  between  the  upper  and  lower  bounds  to  the  sea  brightness  temperature  — 


developed  originally  In  the  context  of  the  geometrical  optics  model  —  is  shown 
to  be  preserved.  The  extensive  numerical  calculations  of  sea  brightness  tempera 
ture  which  were  carried  out  for  the  composite  surface  model  are  summarized,  with 
particular  emphasis  placed  on  the  dependence  on  polarization,  on  frequencies  in 
the  range  1.4  to  20  gHz,  on  observation  angle  between  0°  and  90°  nadir  angle, 
and  on  wind  speeds  from  0  to  20  m/s.  In  addition,  comparisons  are  made  with 
results  based  on  the  single-scatter  geometrical  optics  theory  of  Stogryn,  the 
shadow-  and  double-scatter  corrected  theory  of  Wagner  and  Lynch,  and  the  ex¬ 
periments  of  Hoi  linger. 

In  this  report  we  will  lean  heavily  on  concepts,  equations,  and  results 
developed  in  our  earlier  reports  and  will  refer  to  them  freely  without  re¬ 
definition  or  repetition.  For  example,  the  definition  of  both  the  "Peake 
Integral  method,"  leading  to  the  upper  bound  on  the  t*ve  brightness  temperature, 
and  the  "direct  emission  method,"  leading  to  the  lower  bound,  may  be  found  in 
references  4  and  5;  numerical  results  and  conclusions  based  on  both  methods 
will  be  found  in  reference  6.  In  the  present  study,  both  methods  were  ex¬ 
tended  to  include  small  scale  structure  but  because  of  the  greatt  uracy  of 
the  direct  emission  method  (a  result  established  in  ref.  6)  all  sea  brightness 
temperature  values  quoted  herein  were  calculated  using  this  method.  The  geo¬ 
metrical  optics  equations  also  form  the  basic  mathematical  framework  for  the 
composite  surface  model.  In  this  report  we  develop  the  modifications  to  these 
equations,  but  the  equations  themselves  are  quite  lengthy  and  are  not  repeated 
here;  they  may  be  found  in  ref.  5. 
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II.  COMPOSITE  SURFACE  MODEL 

For  radiation  wavelengths  in  the  centimeter  region  there  will  be  roughness 
of  the  sea  surface  having  scaies  larger  than,  comparable  to,  and  smaller  than 
the  radiation  wavelength.  The  effect  on  the  emissivity  of  swells  and  large 
gravity  waves  can  be  satisfactorily  treated  by  the  geometrical  optics  theory 
while  small  gravity  and  capillary  waves  require  a  perturbation  analysis,  des¬ 
cribed  in  the  next  Section.  Since  only  large  and  small  scale  roughnesses 
(relative  to  wavelength)  are  amenable  to  analysis,  we  assume  the  total  roughness 
spectrum  to  be  divided  into  two  parts  with  each  part  treated  by  the  appropriate 
method.  The  small  waves  ride  on  top  of  the  large  waves  so  that  the  mean  plane 
of  the  random  height  variations  over  a  small  patch  of  emitting  surface  is  not 
the  horizontal  (z  *  0)  plane  but  is  the  randomly  oriented  local  tangent  plane 
of  the  underlying  large  wave.  The  mean  radiation  emitted  into  a  particular 
direction  is  therefore  the  superposition  of  the  mean  intensities  emitted  by 
arbitrarily  oriented  elements  of  surface  area  containing  the  small  waves, 
weighted  by  the  probabilities  of  occurrence  of  the  possible  oreintations 
(i.e.,  by  the  two-dimensional  slope  density  of  the  large-scale  structure). 

The  addition  of  intensities  emitted  by  contiguous  elements  is  justified  by 
the  incoherence  of  thermal  radiation  (note  however  that  in  the  scattering  of 
thermal  radiation  some  degree  of  interference  is  possible). 

It  is  evident  that  the  emission  from  the  composite  surface  is  described 

5 

in  precisely  the  same  way  as  in  the  direct  emission  method  in  the  geometrical 
optics  model,  provided  (1)  the  emissivity  of  an  element  of  surface  area  is 
taken  to  be  that  of  a  patch  of  water  containing  small  scale  roughness  instead 
of  the  emissivity  of  a  patch  of  smooth  water,  and  (2)  the  rms  slope  is  that  of 
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the  large-scale  structure  alone  rather  than  that  of  the  entire  Surface.  It  was 

i 

6  . 

shown  in  an  earlier  report  that  virtually  all  of  the  radiation  reaching  an 
observer  arrives  directly  from  the  surface  with  at  most  a  few  percent  reflected  1 

i  » 

by  the  surface  between  the  time  of  emission  and  observation.  The  small  emission/ 
reflection  contribution  to  the  emissivity  will  be  calculated  for  the  composite  t 
surface  model  as  well.  Since  this  second-order  contribution  is 'a  Small  cor- 

i 

rection  to  the  first  order  ("direct  emission")  term,  we.  may  make  a  convenient  , 
simplification  in  the  calculation  of  the  scattering  by  the  smal 1-scale  structure  % 

of  the  radiation  emitted  by  o;her  portions  o$  the  surface.  Because  the  wave-  ' 

! 

length  is  long  compared  to  the  roughness  of  the  small  component,  the  scattering 

;  ‘  t  1 

by  an  element  of  surface  will  be  sharply  peaked  in  tlje  local  specular  direction 

(l.e.,  relative  to  the  mean  plane  of  the  surface  element)'.  Since  the  nonspecular 

1 

t  • 

component  of  the  scattered  radiation  is  quite  small  we  may  make  the  approximation 

i  i ' 

that  all  the  energy  scattered  by  an  element  of  surface  area  containing  small- 

S 

scale  roughness  appears  in  the  local  specujar  direction.  In  other  words,  we 
ascribe  to  the  surface  element  an  effective  ref lectivi ty*  r^  (i  ■'  the  .polariza¬ 
tion  index)  defined  by 


where  is  the  emissivity  of  the  surface  element  and  contains  the  contributions 
of  the  small  scale  structure  (Sect.  III). 

i 

Thus,  once  is  known  the  emissivity  of  the  total  (cpmposite)  surface 
may  be  accurately  calculated  from  the  geometrical  optics  equations  for  the 
direct  emission  method,  provided  one  makes  the  following  transcriptions.  If 
e*,  r*  denote  the  emissivity  and  reflectivity  of  smooth  sea  water,  $hen,  wherever 
these  quantities  appear  in  the  geometrical  optics  equations,  one  simply  makes 

i  i  ■ 
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i 


t-  i 


£ 

if 


'  A.  ,* 


the  substitutions 


I 


-*■  r 


(2) 


In  addition,  the  rms  slope  of  the  total  surface  is  to  be  replaced  by  the  rms 
slope  of  the  large-scale  wave  structure  as  the  pertinent  roughness  parameter. 


A  slight  revision  in  the  equations  is  also  required  because  the  quantities 
£s’  rs  now  £*ePent^  on  t*ie  direction  of  the  wind,  relative  to  the  directions 

i  I  | 

propagation  and  the  local  surface  normal,  whereas  e  ,  r  of  course  do  not. 


To  the  apparent  emission  temperature,  calculated  in  this  way,  must  be 
added  the  scattered  sky  radiation  to  obtain  the  total  sea  brightness  tempera¬ 
ture.  The  small  scale  surface  structure  will  also  affect  the  angular  Intensity 
distribution  of  the  scattered  atmospheric  radiation.  That  Is,  although  the 
radiation  scattered  by  an  element  of  surface  area  containing  small-scale 
structure  will  be  sharply  peaked  in  the  local  specular  direction,  a  small  non- 
specular  component  will  appear  In  other  directions  as  well.  It  would  be  de¬ 
sirable  to  include  the  diffuse  component  in  the  single  scatter  term,  although 
it  could  well  be  neglected  in  the  smaller  double-scatter  contribution.  This 
could  in  fact  be  done,  although  practical  difficulties  are  encountered.  The 
radiation  leaving  a  given  surface  element,  into  a  fixed  direction,  can  origi¬ 
nate  from  a  large  angular  sector  of  the  atmosphere  because  of  the  diffuse 
component  of  the  scattering  coefficients.  Thus  a  solid  angle  integration  over 
all  source  directions  is  required  for  every  possible  orientation  of  the  sur¬ 
face  element;  averaging  over  all  orientations  then  demands  an  additional 
double  Integration.  Thus  a  4-fold  Integration  (not  counting  those  required 
In  the  evaluation  of  the  scattering  coefficients  themselves)  must  be  performed 
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„,  ■>  gyring*" 


simply  to  calculate  the  single-scattered  atmospheric  radiation.  Doing  so 
would  expand  the  present  program  beyond  its  intended  scope,  so  a  simplifying 
assumption  will  be  made. 

We  assume,  again,  that  all  the  energy  scattered  by  the  small  structure 
on  a  surface  element  is  concentrated  in  the  specular  peak  with  an  intensity 
characterized  by  the  modified  reflection  coefficient  defined  in  Eq.  (1).  In 
this  way,  and  wi th  the  use  of  Eq.  (2),  the  approximate  effect  of  the  small- 
scale  structure  on  both  the  single  and  double-scattered  sky  radiation  may  be 
calculated  using  the  existing  geometrical  optics  computer  programs,  provided 
only  that  is  known  and  that  the  mean  square  slope  of  the  large-scale  struc¬ 

ture  is  used  as  the  geometrical  optics  roughness  parameter.  Some  error  is 
probably  incurred  as  a  result  of  this  approximation.  However,  the  longer  wave¬ 
lengths  (£  a  few  centimeters)  are  tlv  region  of  our  principal  interest  and, 
presumably,  also  the  region  of  maximum  effects  of  the  small-scale  structure; 
fortunately,  it  is  also  a  region  in  which  the  intensity  of  the  atmospheric 
microwave  radiation  is  relatively  low  so  that  the  error  in  the  sea  brightness 
temperature  introduced  by  our  approximation  to  the  scattering  coefficients  is 
also  diminished. 

The  most  accurate  method  of  calculating  the  apparent  emission  tempera¬ 
ture,  and  therefore  also  the  brightness  temperature,  was  shown**  for  the  geo¬ 
metrical  optics  case  to  be  the  "direct  emission  method"  (which  led  to  the 
lower  bound  on  the  true  brightness  temperature).  This  method  will  be  even 
more  precise  for  the  composite  surface  model,  for,  as  we  have  already  shown, 
no  approximations  to  the  scattering  coefficients  are  required  to  calculate  the 
leading  (direct  emission)  term,  while  the  second-order  (emission/reflection) 
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term  is  always  sufficiently  small  so  that  our  approximation  to  the  diffuse 
scattering  by  the  small-scale  structure,  required  in  this  term,  wi i 1  introduce 
a  negligible  error  in  the  emission  temperature.  In  contrast,  the  "Peake  inte¬ 
gral  method,"  which  led  to  the  upper  bound  on  the  true  brightness  temperature 
in  the  geometrical  optics  case,  is  intrinsically  less  accurate  for  several 
reasons.  First  of  all,  all  terms  in  the  expansion  of  the  emissivity  in  terms 
of  increasing  orders  of  multiple  scatter  are  larger  than  the  corresponding 
terms  in  the  same  expansion  of  the  direct  emission  representation;  since  third 
and  higher  orders  of  multiple  scatter  cannot  be  readily  calculated,  their 
neglect  in  both  methods  leads  to  a  significantly  larger  error  when  the  Peake 
representation  is  used  than  for  the  direct  emission  method.  Secondly,  the 
inclusion  of  small-scale  structure  is  more  difficult  to  do  accurately,  for  the 
largest  (single  scatter)  contribution  requires  explicit  evaluation  of  the  com¬ 
posite  surface  scattering  coefficients;  because  of  the  magnitude  of  this  term 
a  simplified  treatment  of  the  diffuse  scattering  [as  in,  for  example,  Eqs.  (1)- 
(2)1  is  not  advisable  and  a  much  more  complex  analysis  is  required.  This  is 
also  true  of  the  second-order  (double  reflection)  term  which  is  an  order  of 
magnitude  larger  than  the  corresponding  (emission/reflection)  term  in  the 
direct  emission  representation.  Finally,  it  must  be  noted  that  since  the 
emissivity  in  the  Peake  method  is  derived  from  an  integration  over  plane  wave 
surface  scattering  coefficients,  the  possibility  of  coherent  interaction  of  the 
fields  scattered  by  the  small  structure  lying  on  different  surface  elements 
of  the  large-scale  structure  must  be  considered;  in  the  direct  emission  method 
this  question  does  not  arise. 

For  these  reasons,  the  composite  surface  model  based  on  the  direct  emis¬ 
sion  method  is  the  one  used  for  the  calculation  of  the  sea  surface  emissivity. 
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However,  for  purposes  of  comparison,  It  was  also  calculated  using  the  Peake 
representation  and  the  approximations  of  Eqs.  (1)  and  (2),  applied  without 
further  justification.  The  method  was  also  used  for  a  test  of  energy  conser¬ 
vation,  assuming  a  perfectly  conducting  surface. 

We  turn  now  to  a  consideration  of  the  sea  surface  spectrum  and  a  more 
precise  definition  of  the  "small"  and  "large"-scale  surface  roughness.  If  the 
height  of  the  surface  relative  to  its  mean  plane  is  given  by  a  random  function 
z  *  5(x,y),  then  the  covariance  function  is  the  ensemble  average 

C(x,y)  =  (?(xo,yo)  S(xQ  +  x,  yQ  +  y))  (3) 

and  the  vector  wave  height  spectrum  is  its  Fourier  transform  which  may  be 
written  in  the  form 


S(K  ,K  )  »  -s-  | dy  dx  C(x,y)cos(K  x  +  K  y)  . 

x  y  2  ,  x  y 

J  J  o 


A  polar  representation  in  K-space  is  convenient;  with  the  transformations 


*  K  cos4> 
Ky  =  K  sin4> 


S(K  .KjdKdK  -  S '  (K,$) KdKd$  .  (5) 

Ay  a.  y 

A  one-dimensional  spectrum  S(K)  and  an  angle-dependent  factor  F(K,$)  may  be 
defined  as  follows,  assuming  the  x  axis  lies  parallel  to  the  wind  direction, 


KS’(M)  =  $(M) 


S(K)F(K,$),  |«|  <  tt/2 

0  ,  |4>|  >  tt/2  . 
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The  anisotropy  factor  is  assumed  to  be  normalized  such  that 


•tt/2 

F(K,*)d*  -  1  .  (7) 

■  -tt/2 

The  composite  surface  is  defined  as  follows.  All  harmonic  components  of 
the  sea  surface  having  wavelengths  less  than  a  few  radiation  wavelengths  con¬ 
stitute  the  small-scale  roughness,  while  all  harmonic  components  with  longer 
wavelengths  are  defined  to  be  of  large  scale.  Thus,  i f  k  *=  2ir/X  Is  the  radia¬ 
tion  wavenumber  and  a  a  parameter  to  be  specified,  then  the  large-scale  spectrum 
is 


(  S(K)  ,  K  <  ctk 
S,(K)  =  <! 

Lo  ,  K  >  ak  (8) 

and  the  small-scale  spectrum  is 


( 0  ,  K  <  ak 

Ss(K)  -  \ 

\  S(K)  ,  K  >  ak  .  (9) 


The  union  of  the  two  is  then  the  total  spectrum 

S(K)  «*  S£(K)  +  Ss(K)  .  (10) 


Two  derived  quantities  are  required.  One  is  the  rms  height,  0,  of  the  small 

scale  structure  which  is  obtained  from 

*00 

a2  -  S(K)dK  ,  (11) 

■  ak 


and  the  other  is  the  rms  slope  In  the  x  and  y  directions  E  ,  E  of  the  large 

x  y 

scale  structure,  given  by 


■ak 


v/2 

E2  -  dK  K2S^(K)  cos2$F(M)  d4> 

'o  -ir/2 

(•ak  fir/ 2 

E2  -  dK  K2S0(K)  sin24>F(K,$)d$  .  (12) 

Jo  J-n/2 

We  note  that  for  optical  frequencies  (k-«°)  Eqs.  (12)  become  the  mean  square 
slopes  of  the  entire  surface  and  their  values  should  agree,  at  a  given  wind 
speed,  v/ith  the  experimental  values  of  the  upwind  and  crosswind  slopes  given 
by  Cox  and  Munk ,  "  ’ 1 2 

The  equilibrium  spectrum  is  fairly  well  known  for  low  wavenumbers  but 
for  larger  wavenumbers,  particularly  in  the  short  capillary  wave  region,  there 
is  still  considerable  uncertainty  in  both  the  wave  number  and  wind  speed  de¬ 
pendence.  Choice  of  a  suitable  spectrum  therefore  involves  some  measure  of 
arbitrariness.  One  form  of  the  vector  spectrum,  continuous  for  all  wave  numbers, 
was  suggested  by  Pierson.  We  found,  however,  that  the  evaluation  of  Eq.  (12) 
for  optical  frequencies  gave  slope  values  which  were  not  in  satisfactory  agree¬ 
ment  with  the  Cox  and  Munk  results.  A  different  form  for  the  one-dimensional 
spectrum,  S ( K) ,  was  recently  suggested  by  Wu,^  who  reanalyzed  the  Cox  and  Munk 
data  in  terms  of  a  logarithmic  dependence  of  mean-square  slope  on  wind  speed, 
and  established  new  spectral  constants  for  a  Phillips-type  spectrum.  The  Wu 
spectrum  has  wind-speed  dependent  cutoffs  at  both  low  and  high  wave  numbers, 
plus  a  discontinuity  at  the  start  of  the  capillary  range  (Kj.^360m  ^).  One  of 

the  three  discontinuities  was  removed  by  adjoining  to  the  low  end  of  the  Wu 
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spectrum  the  form  suggested  by  Pierson.  In  addition,  the  direction  :t 1 ty 
factor  F(K,$)  of  Pierson  was  used  to  form  the  two-dimensional  spectrum  S(K), 
which  is  then 
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■trjSr-T^JTJ- 355=5^2 


S1(K)  2  4  2 

0.08  cos  4*  +  2.56  cos  4>  +  (0.5  +  0.92  cos  4> 


ttK 


2.56  cos\)  exp(-  Tg  g2/K2W**)  ,  0  <  K  < 


S(K)  = 


\  S  (K) 

'  —  (2  +  cos  *} 


KC  <  K  <  Kv 


Ky  <  K 


(13) 


-2 


where  g  =  9*8  msec  ,  W  is  the  wind  speed  in  meters  per  second, 
and 


1 .o7x10_2KcW2, 


S 1 (K)  =  4.6x10_3k"3  exp(-.74g2/w\2) 

S2(K)  *  3.15x10’V3 


(1*0 


These  equations,  together  with  Eqs.  (8)  and  (9)  define  the  "small-scale"  and 
"large-scale"  components  of  the  sea  surface. 

The  parameter  a  may  be  retained  as  an  adjustable  parameter  to  secure  a 
best  fit  to  the  brightness  temperature  data.  We  chose  instead  to  impose  an 
additional  constraint  to  fix  a  and  then  retain  this  value  for  all  the  calcula¬ 
tions.  The  argument  used  to  establish  this  constraint  is  as  follows.  The 
emlssivity  of  a  slightly  rough  surface  can  be  evaluated  through  second  order 
in  the  perturbation  parameter  (ko) .  To  ;he  extent  that  k a  is  sufficiently 
small  that  higher  order  contributions  are  negligible,  this  result  will  be  exact. 
We  therefore  choose  to  include  as  much  of  the  surface  roughness  as  possible  In 
the  small-scale  spectrum  (i.e.,  a  is  to  be  as  small  as  possible)  while  at  the 
same  time  satisfying  the  requirement  that  (ko)  «  1.  Evaluation  of  Eq.  (11) 
for  a  range  of  values  of  a  led  to  the  choice  a  *  0.2.  For  this  a  we  found 


1i* 


2 

(K a)  'v  0.06,  a  value  which  remained  nearly  constant  [because  of  the  predomi- 
nantly  K  3  dependence  of  $(K)]  over  the  range  of  frequencies  and  wind  velocities 
of  Interest  to  us. 

III.  EMI  SSI  VI TY  OF  THE  SMALL-SCALE  ROUGHNESS 

The  solution  of  Maxwell's  equations  for  the  problem  of  a  plane  wave  inci- 
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dent  on  a  randomly  rough  Interface  between  two  media  was  obtained  by  Rice  by 
a  perturbation  method,  for  the  case  In  which  the  wavelength  is  long  compared  to 
the  rms  excursions  of  the  interface  from  Its  mean  (planar)  value.  An  harmonic 
representation  was  assumed  for  the  components  of  the  electric  and  magnetic 
fields  and  for  the  random  boundary  surface.  Application  of  the  differential  equa¬ 
tions  and  the  boundary  conditions  then  led  to  algebraic  equations  for  the 
Fourier  transforms  of  the  fields  on  the  boundary  surface  which  could  be  solved 
for  each  order  of  perturbation.  Explicit  first-order  solutions,  for  both 
horizontal  and  vertical  polarizations  of  the  incident  field,  and  the  second- 
order  solution  for  horizontal  polarization  were  given  by  Rice;  the  second-order 
vertical  polarization  result  was  later  given  by  Valenzuela.^  The  spatial  fields 
can  be  reconstructed  from  the  Fourier  coefficients  and  the  mean  intensity  of 
the  scattered  radiation  evaluated  in  terms  of  integrals  over  the  surface  height 
spectrum.  A  half-space  integration  over  all  scattering  angles,  for  a  fixed 
incident  wave  vector,  then  gives  the  mean  total  energy  scattered  by  the  surface: 
applying  a  statement  of  energy  balance  then  yields  the  mean  surface  emissivity. 

Such  a  procedure  for  calculating  the  emissivity  is  exceedingly  cumber¬ 
some.  However,  the  calculation  may  be  considerably  simplified.  It  is  shown 
in  the  Appendix  that  the  mean  polarized  emissivity,  e^,  may  be  expressed 
directly  in  terms  of  the  Fourier  coefficients  of  the  scattered  field  and  the 
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surface  height  spectrum.  The  equations  derived  there  are  for  the  case  In  which 
the  z  axis  is  normal  to  the  mean  plane  of  the  surface  and  the  observation  direction 
lies  in  the  x-z  plane;  the  horizontal  and  vertical  polarization  components  of 
the  emissivity  have  the  usual  definition  relative  to  this  choice  of  coordinate 
axes . 

It  is  necessary  to  transform  the  expressions  for  the  emissivity  components 

to  an  arbitrary  orientation  of  coordinate  axes,  for  the  "mean  plane"  of  the 

small-scale  structure  is  to  be  an  (arbitrarily  oriented)  element  of  surface  area 

of  the  large-scale  structure.  A  coordinate  system  has  already  been  defined  by 
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the  geometrical  optics  equations:  the  £  dl rection  is  normal  to  the  mean  plane 

of  the  entire  sea  surface,  while  the  £  direction  was  assumed  to  be  parallel  to 
one  of  ihe  principal  axes  of  the  slope  distribution,  that  is,  parallel  to  the 
wind  velocity  vector  W.  Now  a  local  x' ,  y',  z1  coordinate  system  on  a  surface 

A 

element  is  defined  such  rhat  z'  is  the  direction  of  the  normal,  fi,  to  the  ele¬ 
ment  and  y'  is  normal  to  the  plane  defined  by  z'  and  the  wave  vector  k  pointing 
toward  the  observer: 

A  /\ 

z'  ■  n 

y'  »  z'  x  k/| z*  x  k| 

x'  ■  y'  x  z'  .  (15) 

Now,  the  wind  velocity  at  the  surface  is  presumably  parallel  to  the  plane  of 
the  surface  element,  therefore  we  take  the  direction  of  the  wind  velocity  vector 

A 

W,  in  the  local  coordinate  system,  to  be 

W  -  (z1  x  W)  x  z ' / 1  z *  x  W|  .  (16) 

The  angle  between  the  local  wind  vector  and  the  local  x  direction  Is  given  by 


16 


cosij)  «  W  •  x 


A  A 

U  *  •  v  * 


Since  i*  0,  In  general,  the  directional  wave  spectrum  of  the  small  struc 
ture,  given  by  S(K,<J>)  In  Eqs.  (13),  transforms  on  the  surface  element  to 


S^(M)  -  S(M-x)  .  (18] 

If  s  ,s  are  the  slopes  of  the  surface  element  in  the  x  and  y  directions,  k 
x  y  x 

A*  A 

the  x-component  of  k,  and  0  the  angle  between  k  and  the  local  surface  normal, 
then  It  may  be  shown  from  the  preceding  equations  that  the  random  variable  \J> 


is  given  explicitly  by 


_  /  2  2 

s  cos©  /1+s  +S 

kx  +  - -  —r"  -y-  ■■ 

/l+s^+s^  .  /l+x*  sin© 

x  y  y 


Since  the  emissivity  equations  in  the  Appendix  refer  to  a  local  coordinate 
system,  the  spectrum  (denoted  by  W(K)  in  the  Appendix)  must  be  taken  to  be  S^, 
given  by  Eq,  (18) , 

Finally,  the  "horizontal"  and  "vertical"  polarization  components  of  the 
emissivity  also  refer  to  the  local  coordinate  system.  The  transformation  to 
the  fixed  xyz  coordinate  system  of  the  observer  is  also  required  and  is  given 


e  =  a 


2  '->v  j.  h2 

.  e'  +  bj  e' 


v  ,2  v  2  ,h 
e  =  b^e  +  a|  e  , 


where  a^  and  bj  were  derived  in  an  earlier  report  [Eqs. (5. 1 3) - ( 5 - 16)  of  ref.  5] 

This  completes  the  specification  of  the  equations  required  for  the  com¬ 
putation  of  the  emissivity  and  brightness  temperature  of  the  sea  surface  In  the 
composite  surface  model.  Results  of  the  computations  are  discussed  in  the  next 
Section . 
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Before  leaving  the  topic  of  the  small-scale  emisslvity,  we  shall  descr'oe 
some  of  the  numerical  results  for  the  case  of  a  surface  having  only  small  rough¬ 
ness.  In  this  case  the  mean  plane  of  the  surface  is  z  =  0,  the  local  and  ob¬ 
server's  coordinate  systems  coincide,  and  the  emissivity  equations  of  the 
Appendix  apply  without  transformation,  The  emitting  medium  is  assumed  to  be 
sea  water  and  the  spectrum  to  be  the  small-scale  spectrum  defined  earlier, 
with  a  low-frequency  cut-off  of  0.2k.  The  equations  were  evaluated  numerically 
for  observation  angles  of  0-90°  measured  from  nadir,  azimuthal  angles  of  0-90° 
relative  to  the  wind  direction,  wind  speeds  of  up  to  30  m/s,  a  selection  of 
frequencies  between  1.4  and  100  gHz,  and  both  polarizations.  The  parameter  a 
was  fixed,  as  described  earlier,  at  a  value  of  0.2.  The  perturbation  parameter 
(ka)  was  relatively  constant  over  most  of  the  range  of  the  parameters  although 
at  the  higher  frequencies  and  for  large  wind  speeds  it  increased  in  value  to 
a  point  where  the  validity  of  second-order  perturbation  theory  may  be  question¬ 
able  (e.g.,  at  100  gHz,  W  =  14  m/s,  ko  'v  0.6). 

Since  the  flat  surface  emissivity  of  sea  water  is  well  known  we  show  only 
the  roughness  effect,  Ac  =  r.  (rough)  -  e  (flat),  in  the  accompanying  figures. 
Figure  1  shows  the  dependence  of  Ae  on  observation  angle  for  H  polarization, 
v  -  1.4  gHz,  4)  =  0°  (upwind  observation),  and  winds  of  3»  4,  8  and  i4  m/s. 

The  maximum  roughness  effect  shown  would  correspond  to  a  change  in  effective 
emission  temperature  of  about  2eK.  Figure  2  shows  the  same  cases  but  for  cross- 
wind  observation,  *  90°;  a  somewhat  greater  roughness  effect  is  indicated. 

It  was  established  that  for  all  frequencies,  polarizations,  wind  speeds,  and 
observation  angles  the  dependence  on  fetch  angle  is  given  by 

Ac (\p)  ■  Ae(0)cos^  +  Ae (90)  s i n2ip  . 
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Roughness  dependence  of  the  emissivity  of  a  slightly 
rough  surface,  v  =  1.1*1  gHz,  horizontal  polarization, 
upwind  observation. 


I 


Unless  otherwise  Indicated,  the  remaining  figures  will  pertain  to  tfie  case 
iJj  *  0°.  Figures  3  and  4  show  the  roughness  effect  for  H  polarization  at  fre-  , 

quencles  of  8.34  and  .14  gHz,  and  Fig.  5  that  for  100  gHz.,  For  100  gHz,  as  stated 
earlier,  the  perturbation  results  may  not  be  very  accurate  for  W  si  14  m/s;  never¬ 
theless,  It  Is  worth  noting  that  the  curves,  at  least  to  the!left  of  the  peiak, 

i 

are  essentially  the  same  as  those  given  by  the  geometrical  optics  theory  for  the; 
same  surface.  The  dependence  on  frequency  is  shown  In  another  way  In  Fig.  6  for 

I  1 

a  fixed  observation  angle  of  0  ■  55°.  i 

A  more  complicated  behavior  Is  demonstrated' by  the  vertically  polarized 
radiation,  as  is  shown  in  Fig.  7  for  a  frequency  of  ^.34  gHz  and  winds  of  4,  8, 

l 

and  14  m/s.  At  large  nadir  angles  the  rough  surface  would  appear  to  be  much 

!  '  . 

colder  than  the  smooth  surface.  The  roughness  effects  tend  to  increase  with 

increasing  frequency  for  vertical  polarization  also;  in  Fig.  8  i^s  shown,  on  a 

different  scale,  the  region  0  ^  50°  for  a  frequency  of  '19. 35  gHz.  As  in  the 

1 

geometrical  optics  case  for  vertical  polarization  there  is  a  tendency  for  the 
emlsslvlty  to  be  insensitive  to  surface  roughness  for  angles  in  tlpe  neighbor- 

i  ’  1 

hood  of  60°. 

1 

For  the  shorter  wavelengths  the  structure  of  the  spectrum  i tself  ‘begins 

to  have  a  more  prominent  effect  on  the  emissivity.  1  The  peak  in  the  capiMary 

1  •  '  1 

range  is  approximated  by  a  discontinuity  at  K  «*  so  that  asi  the  wind  speed 
increases  the  high  frequency  cutoff,  K^,  of  the  spectrum  will  lie  first  below, 
then  above,  the  discontinuity;  as  a  result  the  emissivity,  obtained  from  an. 
integral  over  the  spectrum,  will  show  a  discontinuity  in  its  derivative  with 
respect  to  W  when  K  =  Kr,  i.e.,  at  W  ^  8.25’m/s.  The  occurrence  of  this  ef- 

V  v 

I 

feet  may  be  seen,  for  example,  in  Fig.  4  which  shows  a  relatively  large  change  1 

1 

in  Aeu  between  wind  speeds  of  8  and  8.5  m/s.  It  is  also  evident  from  Fig.  9  which 
H  •  , 
shows  the  roughness  effect  as  a  function  of  wind  speed  for  fixed  observation 
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Fig.  3.  Roughness  dependence  of  the  emissivlty  of  a  slightly 
rough  surface,  v  ■  8.36  gHz,  horizontal  polarization 
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Fig.  4.  Roughness  dependence  of  the  emissivity  of  a  slightly 
rough  surface.  V  =  lA  gHz,  horizontal  polarization. 


Frequency  dependence  of  the  emissivity  of  a  slightly 
rough  surface  at  fixed  observation  angle  and  various 
wind  speeds. 
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Fig.  7.  Roughness  dependence  of  the  emissivity  of  a  slightly 
rough  surface.  V  *  8.36  gHz,  vertical  polarization. 


angles  of  8  =  0  and  55°.  at  a  frequency  of  1.4  gHz  and  horizontal  polarization. 

At  very  low  wind  speeds  the  rapid  drop  to  zero  of  the  curves  In  Fig.  9  Is  a 
consequence  of  the  wind-speed  dependent  high  frequency  cutoff,  K^,  diminishing 
to  the  point  where  it  is  smaller  than  the  low-frequency  cutoff,  0.2k,  assumed 
for  the  spectrum;  the  roughness  effect  is  then  identically  zero.  This,  of  course, 
is  a  cirect  result  of  our  assumption  of  the  existence  of  a  low  frequency  cutoff 
for  the  small-scale  spectrum  and  Is  not  a  real  effect;  in  reality,  the  curves 
of  Fig.  9  would  presumably  pass  through  the  origin.  We  note  that  in  the  com¬ 
posite  surface  model  this  effect  will  not  occur  since  the  portion  of  the  total 
spectrum  excluded  from  the  small  scale  spectrum  (values  of  S (K)  for  K  less  than 
the  smaller  of  0.2k,  K^)  are  automatically  included  in  the  large-scale  spectrum. 

Under  certain  conditions  the  discontinuities  in  the  spectrum  will  also 

manifest  themselves  in  the  variation  of  Ae  with  0,  for  fixed  wind  speed.  To 

show  this  we  note  that,  for  a  given  frequency  and  wind  speed,  the  integrations 

over  the  fwo-dimensional  spectrum  in  the  emissivity  equations  (given  in  the 

Appendix)  extend  over  a  fixed  annular  region  of  K  ,K  space.  However,  in  these 

x  y 

equations  the  spectrum  is  not  centered  at  the  origin  but  is  displaced  from  it 
by  an  amount  which  increases  with  0  and  k.  Thus  as  0  increases  from  0°  a 
circle  of  discontinuity  moves  across  the  J<  plane  and  eventually  intersects  the 
integration  annulus  and,  as  0  increases  further,  passes  on  through.  The  emls- 
siyity  will  show  an  effect  of  a  discontinuity  in  the  spectrum  for  that  range 
of  0  for  which  an  intersection  exists  between  the  circle  and  the  integration 
annulus.  This  is  illustrated  In  Fig.  10  which  shows  Ae v  vs.  0  at  a  frequency 
of  19.35  gHz.  The  hump  in  the  curves  at  small  angles  can  be  traced  to  the 
discontinuity  at  K  a  associated  with  the  onset  of  the  capillary  range.  At 
the  low  wind  speed  of  4  m/s  the  discontinuity  having  the  principal  effect  is 


^  J* 


that  at  K^;  It  occurs  at  larger  values  of  0  In  this  instance  and  produces  a 


small  peak  at  ^55° .  For  0  50°  there  is  no  contribution  of  roughness  to  the 


emisslvity  since  the  spectrum,  in  this  particular  case,  will  not  overlap  the 


Integration  annulus.  To  what  extent  these  effects  are  model -dependent  is  not 


known  but  should  be  investigated;  in  any  case,  they  appear  to  be  small  and  will 


produce  a  variation  in  apparent  emission  temperature  of  less  than  1°K. 


IV.  SEA  BRIGHTNESS  TEMPERATURES 


Before  applying  and  using  any  model  it  is  advisable  to  test  it  for  in¬ 


ternal  consistency.  A  necessary  condition  for  the  validity  of  a  scattering  or 

k  5 

emission  theory  is  that  it  conserve  energy.  It  has  been  shown  ’that  an  equi¬ 


valent  condition  is  that  the  emisslvity  of  the  emitting  surface  go  to  zero  as 


the  conductivity  of  the  medium  becomes  infinite.  This  condition  is  always 


satisfied  for  the  direct  emission  method  since  the  emisslvity  of  every  surface 


element  (i.e.,  the  emissivity  of  a  slightly  rough  surface)  may  be  shown  to  be 


zero  for  a  perfect  conductor.  While  the  composite  surface  model  using  the 


Peake  representation  was  not  used  for  emissivity  calculations,  except  for  pur¬ 


poses  of  comparison,  essentially  the  same  model  was  used  to  evaluate  the  scat¬ 


tered  atmospheric  radiation.  Calculating  the  emissivity  of  a  perfect 


conductor  leads  to  the  results  shown  in  Fig.  11  [curves  a.)  and  b.)]  for  a 


frequency  of  1.4  gHz  and  a  wind  speed  of  14  m/s.  The  results  are  quite  good 


with  the  maximum  departure  from  exact  energy  conservation  being  0.015  for  the 


scattering  model.  For  comparative  purposes,  the  result  for  the  purely  geo¬ 


metrical  optics  model  is  shown  in  curve  c.).  (The  corresponding  calculation 

! _ ..  i _ _ _ _  .  . _ . _ ... _  .l _ _ _  — - _ ...  a  .  j _ \  ^ 


using  the  non-energy-conserving  Stogryn  theory  was  reported  earlier.)' 
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Two  other  result*;  concerning  the  relationship  between  the  direct  emission 

method  and  Peake  integral  method  are  of  interest.  First,  as  in  the  geometrical 

optics  model,  the  higher-order  multiple  scatter  contributions  to  the  emissivity 

were  found  to  be  an  order  of  magnitude  larger  (for  H  polarization)  in  the  Peake 

method  than  in  the  direct  emission  method,  implying  a  correspondingly  greater 

sensitivity  to  approximations  which  must  necessarily  be  made  in  the  difficult 

calculations  of  the  higher-order  effects.  Secondly,  it  was  proved  theoretically 

for  the  purely  geometrical  optics  theory  that  the  Peake  method  gave  an  upper 

bound  to  the  effects  of  surface  roughness  on  the  emissivity  while  the  direct 

5 

emission  method  gave  a  lower  bound.  Surprisingly,  the  same  relationship  be¬ 
tween  the  emissivities  calculated  numerically  by  the  two  methods  was  also 
found  to  hold,  without  exception,  in  the  composite  surface  model.  Whether  the 
existence  of  the  bounds  in  the  present  model  is  simply  a  consequency  of  the 
specific  details  ascribed  to  the  model,  or  whether  the  proof  of  the  existence 
of  the  bounds  is  actually  valid  outside  the  domain  of  geometrical  optics,  is 
not  known  at  this  time. 

The  roughness  parameter  in  the  geometrical  optics  theory  is  the  rms  slope 
of  the  entire  sea  surface,  as  given  by  the  Cox  and  Munk  equations,  while  for 
the  large-scale  component  in  the  composite  model  it  is  the  rms  slope  of  the 
large  structure  alone.  From  the  definition  of  the  large-scale  spectrum  it  is 
evident  that  as  the  frequency  decreases  the  rms  slope  will  steadily  decrease 
from  the  Cox  and  Munk  value.  The  dependence  of  the  large-scale  rms  slope  on 
wind  speed  is  shown  in  Fig.  12  for  frequencies  of  »,  100,  19.35,  and  1.4  gHz. 

At  low  wind  speeds  there  is  little  difference  in  the  wind  speed  dependence 
except  for  the  lowest  frequencies.  At  W  «  14  m/s,  the  rms  slope  at  19  $}Hz  is 
about  0.67  of  its  geometrical  optics  value  while  at  1.4  gHz  the  ratio  is  MJ.55* 


* 


If  the  small-scale  roughness  had  no  effect  whatever  on  the  emisslvity,  then 
the  frequency  dependence  of  the  brightness  temperature  would  be  described  by 
Fig.  12.  That  is,  the  microwave  sea  brightness  temperature  would  be  accurately 
described  by  the  purely  geometrical  optics  model  but  with  an  effective  rms 
slope  given,  at  a  particular  frequency  and  wind  speed,  by  Fig.  12.  As  Willi  be 
seen  later,  this  is  approximately  what  does  in  fact  occur,  at  least  at  ;t he 

i 

1 

lower  frequencies  where  the  effect  of  the  small  structure  is  not  large.  At 
19  gHz  the  reduction  in  the  large-scale  slope  (from  the  geometrical  optics 
value)  tends  to  be  compensated  for  by  an  increased  effect  of  the  small-scale 
structure  so  that  the  net  departure  of  the  brightness  temperature  from: the  1 

,  .  l 

purely  geometrical  optics  result  is  not  large. 

l  *  ! 

.  ; 

Sea  brightness  temperatures  were  calculated  numerically  for  6  frequencies 

■  i  , 

between  1.4  and  19.35  gHz  water  temperatures  of  283,  293.  and  298°K,  various  1 

1  •  i 

salinity  values  between  fresh  and  very  saline,  and  atmospheric  humidi ty  values 
corresponding  to  dry,  average,  and  humid  atmospheres.  The  dependence  on  all 

i  ! 

the  parameters  except  frequency,  polarization,,  and,  wind  speed  differ  only  in  ! 
negligible  detail  from  that  given’ by  the  geometrical  opti.cs  model.  Since  these 

i 

6  1 

have  already  been  reported  we  shall  here  be  concerned  only  with  the  de¬ 
pendence  on  those  parameters  for  which  the  brightness  temperature  shows  a 
significant  departure  from  the- geometrical  .optics  results. 

I 

•  1 

All  the  values  quoted  in  the  remainder  of  this  report  rfefer,  unless 
stated  otherwise,  to  the  case  of  an  atmosphere  of  average  humidity,  water 
temperature  of  20°C,  salinity  of  35  ppt,  and  a  plane  of  observation  parallel 
to  the  wind  direction.  Since  it  is  not  possible  to  present,  even  in  summary 
form,  the  results  of  all  the  computations  at  all  .frequencies,  only  selected 
results  at  those  frequencies  for  which  data  is  available  wil.l  be  considered. 
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Figures  13  -  18  show  the  dependence  of  the  sea  brightness  temperature  on 
observation  angle  for  various  wind  speeds,  at  frequencies  of  1.41,  8.36,  and 
19.34  gHz  and  both  polarizations.  Comparing  with  the  corresponding  figures  for 
the  geometrical  optics  model  we  see  that  there  is  no  qualitative  difference 
between  the  behavior  predicted  by  the  composite  surface  and  georr<etrlcal  optics 
models,  except  perhaps  for  an  increased  sensitivity  to  w*nd  speed  at  19.34  gHz 
for  wind  velocities  greater  than  8  m/s,  A  detailed  comparison  between  the  two 
models  Is  shown  in  Figs.  19  *  30  for  the  same  frequencies  and  polarizations 
and  wind  speeds  of  8  and  14  m/s.  Except  possibly  for  large  observation  angles, 
the  composite  surface  model  predicts  brightness  temperatures  which  are  some¬ 
what  larger  than  those  calculated  with  the  geometrical  optics  model;  the  dif¬ 
ference  is  about  1°K  at  the  lower  frequencies  up  to  about  4°K  at  the  highest 
frequency.  It  may  be  shown  that  at  the  low  frequencies  the  reduction  in  the 
rms  slope  (from  the  g.o.  value)  associated  with  the  large-scale  roughness 
component  has  a  negligible  effect  on  the  brightness  temperature  (except  at  very 
1  large  angles)  while  the  small-scale  component  adds  a  small  (generally  positive) 
contribution.  For  vertical  polarization  and  large  angles  there  is  a  competition 
.between  a  decrease  in  the  emissivity  due  to  the  small-scale  structure  (cf. 

Sect.  Ill)  and  an  increase  due  to  the  effect  of  the  reduced  large-scale  slope 
on  both  the  emissivity  and  reflected  atmospheric  radiation.  One  would  expect 
that  If  the  angle  dependence  of  the  diffusely  scattered  sky  radiation  were 
explicitly  taken  into  account  (cf.  Sect.  II),  the  brightness  temperatures  at 
nadir  would  increase  further,  due  to  the  large-angle  scattering  of  radiation 
from  the  hotter  regions  of  the  sky,  while  those  near  the  horizon  would  be 
reduced;  this  effect  is  not  expected  to  be  large,  however. 
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F!g.  13.  Horizontally  polarized  sea  brightness  temperatures 
at  1.41  gHz. 


Fig.  14.  Vertically  polarized  sea  brightness  temperatures 
at  1.41  aHz. 


yF%i 


’  *■  ff&irpvTj.’s 


'■v  ■STi't  ^  ,T175,T’ 


0 


Ftg.  18.  Vertically  polarized  sea  brightness  temperatures 
at  19.34  9Hz. 
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Comparison  of  composite  surface  and  geometrical 
optics  brightness  temperatures. 
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21.  Comparison  of  composite  surface  and  geometrical 
optics  brightness  temperatures. 


Fig.  22.  Comparison  of  composite  surface  and  geometrical 
optics  brightness  temperatures. 
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Fig.  23.  Comparison  of  composite  surface  and  geometrical 
optics  brightness  temperatures. 
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Fig.  2k.  Comparison  of  composite  surface  and  geometrical 
optics  brightness  temperatures. 
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Fig.  25.  Comparison  of1 composite  surface  and  gepmetrical 
optics' brightness  temperatures. 
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Fig.  26.  Comparison  of  composite  surface  and  geometrical 
optics  brightness  temperatures. 
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Fig.  28.  Comparison  of  composite  surface  and  geometrical 
optics  brightness  temperatures. 
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Fig.  29.  Comparison  of  composite  surface  and  geometrical 
optics  brightness  temperatures. 
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Fig.  30.  Comparison  of  composite  surface  and  geometrical 
optics  brightness  temperatures. 


Figure  31  shows  the  wind  speed  dependence  of  the  horizontally  polarized 
sea  brightness  temperature  at  a  frequency  of  19.34  gHz  and  an  observation  angle 
of  55°;  AT  In  this  figure  is  defined  to  be  the  difference  between  the  brightness 
temperature  at  wind  speed  W  and  that  at  zero  wind  speed.  The  solid  curve  Is 
that  calculated  with  the  present  theory,  while  the  experimental  points  were 

g 

reported  by  Hoi  linger.  The  error  bars  on  the  right  of  figure  represent  the 
expected  uncertainty  in  the  data,  as  quoted  by  Hoi  linger.  Experiment  and  theory 
appear  to  agree,  within  the  experimental  error. 

In  Fig.  32  is  shown  the  change  in  the  percentage  polarization,  P,  with 
wind  speed,  where  P  =  100  x  (Ty"T^)/(Ty+T^) .  The  solid  line  represents  the 
theoretical  values  while  the  dashed  line  represents  Hollinger's  linear  fit  to 
his  data.  The  two  appear  to  be  in  agreement  over  the  range  of  mearured  values. 
However,  if  the  percentage  polarization  P,  rather  than  AP,  is  compared  with 
experiment,  a  constant  difference  of  several  percentage  points  is  revealed, 
even  for  zero  wind  speed.  Although  a  slight  swell  may  remain  at  zero  wind 
speed,  a  simple  flat  surface  calculation  —  independent  of  any  theory  or  model 
of  surface  roughness  effects  —  should  nevertheless  correspond  quite  closely  to 
the  measurements.  This  is  not  the  case  and  the  reason  for  the  discrepancy  is 
not  known,  although  a  slight  polarization  mixing  in  the  antenna  might  be 
suspected. 

In  Figure  33  we  show  the  wind  speed  dependence  of  the  vertically  polari¬ 
zed  brightness  temperatures  at  a  fixed  observation  angle  of  55°  and  for  the 
three  frequencies  of  1.41,  8.36,  and  19-34  gHz.  For  the  corresponding  measure¬ 
ments  of  Hoi  linger,  which  are  not  shewn  in  the  figure,  we  refer  the  reader  to 
reference  8.  Hoilinger  fit  his  data  with  straight  lines  and  quoted  values  of 
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19.34  GHz 


Fig.  33.  Wind  speed  dependence  of  the 


the  slopes  of  the  lines  at  the  three  frequencies.  Because  of  the  obviously 
nonlinear  behavior  of  the  theoretical  curves  It  is  difficult  to  assign  a 
single  value  to  d(AT)/dW.  However,  examination  of  the  data  points  and  their 
expected  uncertainty  show  the  theoretical  curves  o.:  Fig.  33  to  lie  well  within 
the  range  of  the  measurements. 


V. 


SUMMARY  AND  CONCLUSIONS 


The  second-order  geometrical  optics  theory  of  emlsslvlty  and  brightness 
temperature  of  a  rough  surface  is  extended  to  low  frequencies  by  the  Intro¬ 
duction  of  a  composite  surface  model  In  which  the  diffraction  effects  of  surface 
roughness  smaller  than  a  wavelength  are  explicitly  included.  The  total  surface 
height  spectrum  is  divided  into  a  large-scale  and  a  small-scale  spectrum,  with 
the  separation  point  determined  by  the  radiation  frequency.  The  small-scale 
spectrum  is  assumed  to  describe  the  surface  structure  residing  on  a  surface 
"element'1  of  the  large-scale  structure.  Equations  for  the  emlsslvlty  of  the 
small  structure  are  derived  using  second-order  perturbation  theory.  This 
roughness-dependent  emisslvlty  replaces  the  elemental  smooth  surface  emlsslvlty 
in  the  geometrical  optics  equations;  the  effect  of  the  large-scale  roughness 
Is  described  by  the  same  equations  but  with  the  roughness  parameter  now 
determined  by  only  the  large-scale  spectrum.  The  composite  surface  model  Is 
shown  to  conserve  energy  to  a  high  degree  of  approximation.  The  "direct  emis¬ 
sion"  method  of  calculating  rough  surface  emissivity  is  used  to  evaluate 
horizontally  and  vertically  polarized  sea  bright  brightness  temperatures  over 
a  wide  range  of  frequencies,  wind  speeds,  water  temperatures,  salinities,  and 
atmospheric  conditions.  The  results  are  compared  with  those,  of  the  purely 
geometri  al  optics  model  and  the  differences  are  shown  to  be  no  more  than  a 
few  degrees  k°K)  except  for  very  large  observation  angles.  Significant 
differences  18K)  do  occur  for  all  frequencies,  both  polarisations,  and  most 
wind  speeds.  Thus  for  applications  which  require  calculations  of  this  accuracy, 
use  of  the  composite  surface  model  is  indicated.  Theoretical  and  experimental 
sea  brightness  temperatures  for  both  polarizations,  three  frequencies,  and  a 
range  of  wind  speeds,  are  in  good  agreement. 
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APPENDIX  —  EMISSIVITY  OF  SLIGHTLY  ROUGH  SURFACES 

In  this  Appendix  we  derive  the  equations  for  the  emissivity  of  a  rough 
surface  whose  height  variations  are  small  compared  to  the  radiation  wavelength. 
Expressions  for  the  scattered  electromagnetic  field,  defined  as  the  difference 
between  the  total  field  and  that  due  to  a  smooth  plane,  are  used  in  a  statement 
of  energv  conservation  to  yield  an  "optical  theorem"  for  rough  surface  scatter¬ 
ing.  This  equation  expresses  the  emissivity  as  a  sum  of  the  flat-surface 
emissivity  plus  terms  depending  on  roughness.  Use  of  a  Fourier  integral  repre¬ 
sentation  for  the  scattered  fields  then  provides  an  equation  in  which  the  emis¬ 
sivity  is  expressed  directly  in  terms  of  the  Fourier  coefficients  of  the 
scattered  fields.  Perturbation  solutions  for  the  Fourier  coefficients,  known 
through  second  order  in  the  perturbation  parameter,  are  inserted  into  the 
emissivity  equation,  and  averaging  over  an  ensemble  of  surfaces  having  random 
surface  height  variations  then  leads  to  the  desired  equation  in  which  the 
emissivity  is  given  explicitly  in  terms  of  the  two-dimensional  surface  height 
spectrum  of  the  surface. 

I .  Asymptotic  Limit  of  the  Field 

We  consider  a  plane  wave  incident  on  the  rough-surface  boundary, 
z  =  s(x,y) ,  which  separates  the  source  medium  from  a  dielectric  with  index  of 
'•^fraction  n  (Fig.  34).  The  electric  field  satisfies  the  vector  He’mholtz 
equation 

7  x  V  x  E  -  k2E  =  0  (1.1) 

everywhere  above  the  surface.  We  first  consider  the  case  of  horizontal  polari¬ 
zation.  Then,  if  the  direction  of  incidence  is  chosen  to  be  in  the  x-z  plane, 
the  following  Fourier  representation  for  the  electric  field  applies: 

Eh(r)  =  y(exp(iko*r)  +  Rh(cos6o)exp( ikQ*r) ]  +  E^c ( r)  ,  (1.2) 
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where 


!sc  “  dKx  ^Ky  Ah(Kx,Ky)exp[!(Kx  x  +Ky  y+bKz)}  ; 


(1.3) 


-00  4  -00 


bK  -  (k2-K2-K*> 
bK  •  i('<x+'<y-k2),S 


9  9  * 

(K  +K  )  <  k 

x  y 

2  2  ^ 

.(iC+K  )  >  k 

■  x  y  i 


(1-4) 


k  =  k(-sfn8  x  -  cos0  z)  , 

~o  o  o 

k'  a  k(-sln0  x  +  cos0  z) 

~o  o  o 


(1.5) 


The  first  term  in  the  expression  for  the  field  is  the  flat-plane  scattering 
result  with  reflection  coefficient  R^.  A  perturbation  expansion  for  the  Fourier 

coefficients,  valid  when  both  the  surface  slope  and  the  ratio  of  surface  height 

! 

1 5 

to  wavelength  are  small,  was  derived  ,by  Rice,  and  thq  first  few  terms  may  be 

i  , 

found  in  Sec.  1 1 !  of  this  Appendix,  1  ‘ 

.  i 

The  asymptotic  form  of  Eq . , (1.2)  determines  the  scattering  coefficient. 


We  first  introduce  .polar  coordinates  in  Eq.  (1.3)  and  then  write 


where 


E*1  -  J^1  +  J^1 

~sc  XI  +  X2 


h  h 

=  dK  K  du  A  (Kcosco,  Ks.Inu)  tex|b  iK(xcosuH-ysinu))  ]  x 

V°  0  '  !  ,  9  h 

xfexp  i  (k  -K*),  z] 
r00  ir  ... 

=  dK  K  doj  A^1  (KcosUjKs  inw)  [exp  i'K(xcosw+ysinu>)  ]  x  1 


(1.6) 


(1.7) 


x  exp[-(K^7k^)  z]  .  ,  (1.8) 
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We  change  Integration  variables  In  Eq.  (1.7)  via  the  transformation  K  *  kslny 

and,  also,  we  introduce  spherical  coordinates  (r,0,4>)  for  (x,y,z)  : 

1  2tt  /•tt/2 

=  k^  da)  dy  slny  cosy  A^(ksInycosu),ksinysinU))  exp[  1  kr(k*  r)  ]  , 


where 


k*r  =  cosy  cos0  +  siny  sln0  cos(to*'<j)) 


(1.9) 


(1.10) 


The  following  relation  '  between  plane  wave  and  spherical  waves  holds  in  the 
asymptotic  limit: 

exp[  i kr (k* r)  ]  (-—)  [6 ( 1  -k« r)  exp(ikr)-6(1+k*r)exp(-Ikr)  ]  .  (1.11) 

By  substitution  of  Eq.  (1.11)  into  Eq.  (1.9),  we  then  obtain 


J*?  —*■  (~— •)  (cos0)  Ah(ksin0cos<J>,ksln0sin<|>) 

—  i  r-* 00  i  ~  r 


(1.12) 


for  the  delta  function  which  multiplies  the  incoming  spherical  wave  is  always 
zero  for  the  integration  range  0  <  y  <  tt/2 . 

h 

The  remaining  quantity,  J,j»  contains  a  dying  exponential  factor.  Because 

h 

of  this  exponential  factor,  the  magnitude  of  can  be  shown  to  drop  off  at  least 
1  _2 

as  rapidly  as  z  ,  when  z  is  more  than  a  few  wavelengths  above  the  medium  plane. 

1  h 

The  dominant  term  in  the  asymptotic  region  is  therefore  Jj.  Thus,  from  Eqs. 

(1.2),  (1.6),  and  (1.12),  we  have 

'  1^0  rV»  YtexpOko-r)  +  Rh  exp (ik^r)  ]+Fh (k ,  kQ) ,  (1.13) 


where 


p^(k,ko)  ■  (^p)  (cos0)  Ah(ksin0cos<{>,  ksin0sin<f>)  . 


(I.IM 


V  ^ 


where  S  represents  the  hemisphere,  p  is  the  rough  surface,  and  T  is  the  total 
surface.  Equation  (2.1)  is  a  statement  of  energy  conservation;  i.e.,  absorption 
and  scattering  by  the  rough  surface  must  balance  the  incident  energy  flow.  This 
means  that  the  left-hand  side  of  Eq.  (2.3)  must  be  proportional  to  the  incident 
power  minus  the  total  scattered  outflow.  Indeed,  the  substitution  of  the  usual 
asymptotic  form  of  the  scattering  solution, 


r-Tw 


sxp(ik  *r)  +  fh(k,k  ) 


(2.4) 


(2.5) 

(2.6) 

(2.7) 


The  quantity  y  is  the  scattering  coefficient;  it  is  defined  as 


So  Y  <M0> 


f  (k.kQ) r 


(2.8) 


Here,  the  power  incident  on  the  segment  of  rough  surface  subtended  by  the  hemis¬ 
phere  S  (of  radius  r)  is  unit  intensity  times  Sq,  the  projection  of  the  surface 
area  in  the  direction  of  incidence: 

S  ■  irr^  cos0  .  (2.9) 

o  o 

h 

Finally,  en(-k  )  is  the  horizontally  polarized  emissivity  of  the  rough  surface, 

2 

with  Eq.  (2.7)  following  from  Eq.  (2.6)  via  a  theorem  due  to  Peake.  Equation 
(2.7)  can  be  written  explicitly  in  terms  of  the  fields  by  use  of  Eq.  (2.2): 
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eh(_bo) 


-u  C  h 

-g2_  Re  dS*(EhxHh  ) 
o 

;S 


(2.10) 


We  now  proceed  to  the  derivation  of  the  optical  theorem  by  the  substi¬ 
tution  of  the  asymptotic  representation  Eq.  (1.13)  into  Eq.  (2.10).  This  ap- 

U 

proach  Insures  a  nonzero  interference  term.  The  asymptotic  magnetic  field,  H  , 
is  determined  from  E^1  via  Maxwell's  equations: 


yoc  ~  r+a>  ^k0x  y)exp(iko*r)  +  (k^x  y)  Rh  exp(iko*r) 


+  x  Fh(k,k  )]iSi£iMi 


(2. to 


Three  sets  of  terms  then  arise  from  the  cross  multiplication  in  Eq.  (2.10).  The 

first  set,  involving  products  of  plane-wave  terms  only,  must  lead  to  the  flat- 

h 

plane  smissivity  as  roughness  effects  are  contained  only  in  F  .  The  second  set 
of  terms  describes  interference  between  the  specularly  scattered  plane  wave  and 
the  outgoing  spherical  wave,  with  the  result  that  the  scattering  coefficient  is 
evaluated  in  the  specular  direction.  The  third  set  contains  the  single  term 
representing  the  scattered  intensity.  After  some  algebra  and  frequent  ap¬ 
plication  of  Eq.  (1.11),  we  arrive  at  the  following  result: 


'(-k0)  1  -  |  Rh(cos0Q)  |  2  +  (jg-  j-^n{R^(coseo)[yFh(k^,ko)]} 


-  ~  dn|F 

r\  J . . 


Fh(kr,ko)|  2 


(2.12) 


The  angular  integration  is  over  the  upper  hemisphere. 

The  first  term  in  Eq.  (2.12)  is  the  emissivity  of  a  plane  surface.  The 
remaining  two  terms  describe  the  effects  of  surface  roughness.  A  check  on  this 


result  follows  from  consideration  of  a  perfectly  conducting  surface.  Then,  by 
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Kirchhoff's  law,  the  emissivity  must  vanish  because  the  absorptivity  is  2ero. 
Since  R^  =  -1  for  a  perfect  conductor,  Eq.  (2.12)  reduces  to 


JmtyFh(k^,kJ  ] 


~o  ~o 


dn|Fh(kr,kQ)|2 

JU 


(2.13) 


This  is  the  expected  optical  theorem  for  a  conductor;  the  minus  sign  in  Eq . (2 . 13) 
arises  because  it  is  the  reflected  plane  wave  which  interferes  with  the  spheri¬ 
cal  wave. 


Finally,  we  relate  the  emissivity  to  the  Fourier  transform  of  Sec.  I 
by  substitution  of  Eq.  (1.14)  into  Eq.  (2.12).  Because  of  the  choice  of  the 
magnetic  field  to  describe  vertical  polarization,  Eq.  (1.15),  the  same  formal 
relation  between  emissivity  and  Fourier  transform  holds  for  both  polarizations: 

2 

e'(-k)  ■  [1-|  R.  (cos0  )  |  2]  -  (cos0  )Re/Rj  [yA1  (k  ,0) 

"»U  1  v*  O  O  *  I  "**  OX  J 

0 

2 

-  dQ  cos^0|a'  (ksin0cos<}>,  ksin0sin<f>)  |  , 

°  U 

i  *>  h,v  (2.14) 

Rough-surface  corrections  to  the  emissivity  can  now  be  determined  from  approxi¬ 
mate  expressions  for  the  Fourier  transform. 


Ill .  Ensemble  Averages 

We  can  write  the  vector  transform  in  terms  of  components: 

A'(Kx,Ky)  =  aiCVKy)  ^  +  31  (Kx,Ky)  y  +  Y!(Kx,Ky)  z  ,  i  =  h,v.  (3.1) 

In  the  procedure  due  to  Rice,  each  component  Is  expanded  as,  say, 

a'(Kx,Ky)  =  a|(Kx,Ky)  +  a2(Kx,Ky)  +  .  *  I  -  h,v  ,  (3.2) 
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where  <xj  is  first  order  in  surface  height,  is  second  order  in  surface  height, 

etc.  The  a!  are  determined  by  the  requirement  of  continuity  for  the  tangential 

electric  and  magnetic  fields  at  the  surface  boundary.  Explicit  expressions  for 

the  analogous  first-  and  second-order  Fourier  series  coefficients  •/ere  given  by 
1 8  ' 

Valenzuela.  These  expressions  have  been  transformed  here  to  reflect  nots- 
tl-jnal  differences  as  well  as  our  use  of  Fourier  integral  rather  than  Fourier 


series. 


The  results  for  horizontally  polarized  incident  radiation  will  be  con- 

h  h  h 

sidered  first.  The  components  (a,  6  ,  y  )  correspond  directly  to  the  Valen- 
h  h  h 

zuela  (A  ,B,Cn).  If  the  spectrum  of  the  surface  profile  is  represented  as 


C(x,y)  *  dp  dq  S(p,q)  exp[i(Pxx  +  qy)]. 


(3.3) 


S  (p,q)  =  S(-p,-q)  , 


(3.3a) 


then  the  first-order  coefficients  are 


ot,  (K  ,K  )  = 

1  x’  y 


y,  S(K  -k  ,K  ) 
1  x  ox  y 


8,(K  ,K  )  =  y"  S(K  -k  ,K  ) 

1 '  x’  y  2  x  ox’  y 


Y?(K  ,k  )  = 

1 1  x*  y 


y,  S(K  -k  , K  ) 
3  x  ox’  y 


(3-M 


where 


■IV.  K  K 
h  x  v 


(c+bn2) 


<\^V 

(c+bn2) 


-iV. CK 

_ h...  1. 

(c+bn2) 


(3.5) 


Vh  «  (nZ-1) (1+Rh)  , 

b  =  (k^K^-K2)55  , 

1^ 

„2  2  „2  ,.2.  2 

c  -  (k  n  -Kx-Ky)  . 


(3.6) 


Inspection  of  Eq.  (2.14)  shows  that  only  B,(K  K  )  is  required  to  provide  an 

^  x  y 

expression  for  emissivity  valid  thru  second  order  in  surface  height.  Indeed, 

we  require  B,  evaluated  at  K  »  k  *»  ksin0  ,  K  =0: 

«•  x  ox  o  y 

»  r  r 

(cos0o.(n2-sln\) ^B^k^.O)  •  dKx  dKy [v^ ( , Ky )]  Stk^-K,,-^) 

•00  -00 

*  s<Vkox-V’  (3-7) 


where 


vb  =  k2Vh(n2-sin20o) 

h  2  (bc+Kx} 

v»  «  -k(n2-l)V. - 5- 

2  h  (c+bn2) 


(3.8) 
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The  approach  for  vertically  polarized  Incident  radiation  must  be  modi¬ 


fied  slightly  because  the  coefficients  (AV,BV,CV)  as  given  by  Valenzuela  are 


for  the  electric  field  while  the  coefficients  (av,0  ,y  )  of  this  paper  are  for 
the  magnetic  field.  However,  It  Is  easy  to  show  that  the  two  sets  of  coeffi¬ 


cients  are  related  by 


“ X’V  - 


eV(Kx,Ky)  =  k’1(bAV-KxCV) 

yv(Kx,Ky)  =  k_1(KxBv-KyAV) 


■(kb)‘,[KxRAv+(bV)Bv]  , 


(3.9) 


where  (A  ,B  ,C  )  are  functions  of  (K  ,K  ).  We  can  now  write  down  the  flrst- 

x  y 

order  coefficients  of  interest: 


“X-V  *  ^  s(vkox-V  • 

e;(Kx,K)  -  V\  S(Kx-kox,K >  , 


yY(K  ,K  )  -  p*  S { K  -k  ,K  )  , 

'  1  x’  y  3  x  ox  y 


(3.10) 


where 


IV  K  K  o  o  4 

M,  =  “koy  +  X  (n2'sln20o)  (c-b) 

1  (c+bn2)-  ox  kn2  0 


iV  (bc+K2+K2)  fK 
v'  x  y'  ' 


(c+bn2) 


(n  -l)(l+Rv) 


fK  ,  2 

(n2-sln20  )  , 

Lkn*  J 


(3.11) 


(3.12) 
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The  second-order  coefficient,  ^^ox*^’  reclu'red  f°r  the  emissivity  calcula¬ 


tion  is 


[n2cos0o+(n2-sin20o)  2]$^{ko;,,0)  *  dKx  dKy  Iv^+V^K^Ky)]  x 

1  -  J  -00 


where 


x  S (k  -K  ,-K  )$(K  -k  ,K  )  ,  (3.13) 

ox  x’  y'  x  ox  y  * 


9  9  9  "2 

=  -k  v  (n  -s I nz0  ) 
v  o 

=••  2k  V  (n2-sin20  ?  *  K 

OXV  0  (c+bn2)  x 


kV  (n2-1)  (n2-sln20 J  .  ,99 

- 0  (bcfK2}  -s,n2e  (k2+x2) 


(c+bn2) 


o  x  y' 


(3.14) 


A  problem  of  considerable  interest  Is  that  for  which  only  statistical 
information  about  the  surface  is  available.  We  will  consider  stationary  random 
surfaces,  so  the  average  height  ^?(x)^  =  0,  and  this  implies 

<A|  >  =  (A*)  =  0  .  (3.15) 


The  first  approximation  to  the  ensemble  average  of  Eq.  (2.14)  is  then 

(ej(-kQ))  -  [l-|K,(cos0o)|2.|-'T"  (cos0o)Re[R'‘(cos0o)(e2(kox.°))] 


2  * 

_  (2irk)__  dQ  cos2q  ^  |a|  (ks i n9c0S(J> , ks I n©s i n(^)  |  ^  , 
0 


I  *  h^v 


(3.16) 


Thus,  the  correction  to  the  flat-surface  emissivity  is  second  order  In  surface 
height. 
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The  expectation  values  for  Eq.  (3.16)  are  easily  obtained.  Because  of 
the  restriction  of  stationarl ty,  the  correlation  function  of  the  random  process 


(?(x,y)  ?(x  +  T.,  y  +  tj)  ■  02C(t)  , 


(3.17) 


2  2 

where  T  *  (ti+t2^  *  can  s^own  6y  use  of  the  inverse  of  Eq.  (3-3)  that 


(s(pj  ,q1)S(p2,q2))  ■=  W(p2,q2)  ~  dp  p  JQ  p/(p1+p2) 2+(q^q2) 2 


(3.18) 


Here,  the  energy  spectrum  W  is  the  Bessel  transform  of  the  correlation  function, 


W(pl»q1)  *  2tt02  dt  t  C(t)  Jo|x(p2+q2)  ]  . 


(3.19) 


The  other  factor  in  Eq.  (3.18)  is  obviously  a  delta  function  but  we  leave  It  in 
integral  form  because  the  required  expectation  values  necessitate  the  choice 
Pj  *  -p2>  q^  =  -q2.  The  infinity  arises  because  of  the  representation  of  the 
incident  beam  by  a  plane  wave.  However,  we  can  relate  Eq.  (3.18)  to  the 
(infinite)  beam  size  Sq  defined  in  Eq.  (2.9)  by  the  intermediate  step 


■r 

^S(-p,-q)S(p,q) )  =  W(p,q)  •—  lim  dp  p 

r->  oo 

•ft 


=  W(p,q) 


lim  r  cos0 


o  r  + 00 


=  W(p,q)  ~4  -=§-  • 

(2tt)  2  cos0o 


(3.20) 


We  can  now  write  down  the  expectation  values  required  in  Eq.  (3.16) 
From  Eq.  (3.*0,  we  have,  say, 
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(3.21) 


(la1?!  )  •=  | y1? | 2  (s,V(K  -k  ,K  )  S(K  -k  ,K  )) 

\  1  '  1  ^  1 1  \  x  ox  y  x  ox*  y  ' 


h  2  s 
IpVI  W(K  -k  ,K  )  - T 

"i 1  x  ox’  y  /.  .2 


(2it)  cos0. 


,  2  h  2 

and  similarly  for  the  averages  of  |0,j  and  |>'j  |  •  The  corresponding  result  for 


vertical  polarization  is 
2\  .  2 
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(3.22) 


The  second-order  averages  again  follow  by  use  of  Eq.  (3.20) 
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dK 


dK  f vY+v^ ( K  ,K  )  w(K  -k  ,K  )  . 
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(3.24) 


The  area  Sq  in  Eqs.  (3.21)  -  (3.24)  cancels  the  like  quantity  in  Eq.  (3.16). 
Examination  of  the  above  results  demonstrates  that  the  rough-surface  effects  on 
emissivity  are  determined  by  the  energy  spectrum  for  the  surface. 
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